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ABSTRACT 


The  use  of  a  linearized  mathematical  spline  for 
Interpolation  between  given  points  occasionally  yields 
extraneous  Inflection  points  which  for  some  applications, 
notably  the  fairing  of  ship  lines,  are  not  acceptable. 

To  avoid  this  difficulty,  a  spline  in  tension  is  con¬ 
sidered  and  the  value  of  tension  which  is  sufficient  to 
remove  all  the  extraneous  inflection  points  is  determined. 
The  slopes  of  the  interpolating  curve  at  each  of  the  ^ 
given  stations  are  considered  as  the  unknowns  and  can  be 
calculated  without  the  necessity  of  numerical  matrix 
inversion. 
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I.  INTRODUCTION 

Frequently  In  engineering  work  one  encounters  the 
need  to  relate  a  smooth  curve  to  a  large  number  of  discrete 
points.  Where  the  discrete  points  represent  experimental  data^ 
one  usually  desires  a  smooth  curve  which  vflll  approximate  the 
given  data;  thus,  the  cxirve  Is  not  required  to  pass  through 
every  given  point.  However,  when  the  discrete  points  are  In¬ 
tended  to  Indicate  the  geometric  outline  of  a  physical  body, 
the  smooth  curve  Is  usually  required  to  pass  through  the  exact 
locations  of  all  the  given  points;  such  a  curve  serves  the 
purpose  of  Interpolating  points  between  the  given  points. 

Only  the  second  type  of  curve,  the  Interpolation  curve,  will 
be  considered  here. 

Historically,  draftsmen  have  determined  a  "smooth", 
or  "faired",  interpolation  curve  with  the  aid  of  a  long  flexible 
beam  or  "spline".  The  spline  is  constrained  to  pass  through 
the  given  points  In  a  plane,  and  heavy  objects  -  "dogs"  -  placed 
between  the  given  points  keep  the  spline  from  moving  while  Its 
deflection  curve  is  being  outlined  on  the  drawing.  Although 
seemingly  not  subject  to  variations,  two  draftsmen  will  seldom 
achieve  Identical  results.  In  spite  of  these  slight  variations 
the  method,  in  general,  yields  highly  saclsfactory  results. 

The  usual  mathematical  analysis  of  the  draftman's 
spline  considers  the  small  deflections  of  a  simply  supported 
elastic  beam.  Using  "Strength  of  Materials"  assumptions,  the 
deflection  curve  can  be  represented  by  cubic  polynomials 
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connecting  the  given  points  with  adjacent  cublcs  having  their 
first  and  second  derivatives  equal  at  their  common  given  point. 
The  method  has  been  aptly  termed  a  "piecewise  polynomial"  method 
by  Blrkhoff  and  Garabedian  [l];  a  derivation  of  the  fom  used 
by  them  for  one  variable  (the  slopes  at  the  given  points  are 
taken  as  the  unknowns)  appears  In  Part  II.  Although  Schoenberg 
notes  (2]  related  work  on  piecewise  polynomial  approximation 
found  in  actuarial  literature  on  "osculatory  interpolation", 
he  Is  apparently  the  first  to  extensively  treat  the  polynomial 
analog  of  the  draftsman's  spline  [3»^J. 

The  cubic  polynomial  "spline"  curve  has  many  useful 
properties  [5,6,7]  and  has  been  successfully  used  as  an  Inter¬ 
polating  curve  In  many  problems  far  removed  from  those  treated 
by  Its  physical  counterpart,  the  draftsman's  spline.  Unfortu¬ 
nately,  the  cubic  polynomial  method  Is  not  satisfactory  In  the 
problem  closest  to  Its  origin:  producing  a  smooth  geometrical 
outline  which  passes  through  the  given  points .  The  difficulty 
centers  around  ceirtaln  conditions  on  curvature  implied  by  the 
physical  problem.  One  may  require  a  convex  interpolation  curve 
for  a  "convex"  set  of  given  points.  For  instance,  consider  the 
five  points  In  Figure  1;  the  Interior  points  2,  3,  and  4  have 
positive  second  differences  (see  eq.  22)  and  it  would  be 
reasonable  to  require  the  Interpolating  curve  to  have  a  positive 
second  derivative.  Such  a  requirement  Is  indispensable  where  the 
Interpolation  curve  represents  the  geometrical  outline  of  a  form 
which  must  be  fabricated,  such  as  a  ship's  hull  or  automobile 
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body.  Hie  cubic  polynomial  method  does  not  always  satisfy  this 
plausible  requirement;  under  conditions  discussed  in  Part  2.2, 
the  curve  will  have  a  negative  second  derivative  in  certain 
regions  where  a  positive  second  derivative  is  expected  from 
the  location  of  the  given  points,  thus  Introducing  unwanted,  or 
extraneous.  Inflection  points. 

Customarily,  the  draftsman’s  results  do  not  have  the 
extraneous  Inflection  points  predicted  by  cubic  polynomial 
representation.  Thus,  by  refining  the  assumptions  made  to  obtain 
the  cubic  polynomial  representation,  we  may  possibly  obtain  a 
better  mathematical  representation  which  duplicates  the  physical 
model's  lack  of  extraneous  inflection  points.  The  assimiptlons 
can  be  considered  as  follows:  l)  the  group  of  "Strength  of 
Materials"  assumptions  concerning  the  distribution  of  stress  and 
strain  in  the  beam;  2)  small  deflections,  or  more  precisely, 
curvature  taken  equal  to  the  second  derivative;  3)  simple 
supports  at  the  given  points .  The  validity  of  the  first  group 
of  assumptions  seems  certain,  at  least  with  respect  to  extraneous 
inflection  points,  since  experiment  or  physical  intuition  can 
quickly  supply  examples  of  simply  supported  physical  splines 
which  have  extraneous  inflection  points  under  approximately  the 
same  conditions  as  those  which  yield  such  points  in  the  cubic 
polynomial  method.  The  second  assumption  can  be  dismissed  from 
consideration  since  it  will  be  shown  that  the  conditions  under 
which  an  extraneous  inflection  point  occurs  in  the  cubic 
polynomial  method  is  independent  of  the  magnitude  of  the 
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deflections.  The  assumption  of  simple  supports,  however,  is 
Imned lately  open  to  question.  The  heavy  "dogs"  provide  the 
physical  mechanism  which  could  produce  reaction  moments  and 
maintain  tension  in  the  spline.  It  seems  likely  that  the  drafts 
man  recognizes  the  "convexity"  of  the  given  points  and  manipu¬ 
lates  the  spline  by  deviating  x’rom  simple  support  conditions  in 
order  to  achieve  the  implied  curvature  conditions,  thereby  avoid 
ing  extraneous  inflection  points. 

The  question  now  arises  as  to  what  parts  of  the  simple 
support  assumption  can  be  successfully  retained.  Since  only  a 
siaall  eimount  of  tension  applied  to  a  simply  supported  spline 
would  physically  remove  the  extraneous  inflection  points  in 
most  cases,  one  is  led  to  consider  the  simpler  condition  of 
tension  alone  and  neglect  the  reaction  moments . 

The  value  of  tension  can  be  expected  to  change  abruptly 
at  each  of  the  "dogs".  However,  as  a  mathematical  expediency, 
tension  will  be  considered  constant  over  the  length  of  the 
spline.  This  assumption  permits  the  use  of  a  mathematical  model 
in  which  the  tension  is  produced  by  axial  loads  at  the  end  of 
the  spline  rather  than  by  reaction  forces  at  the  support  points. 
Thus,  the  net  result  of  these  changes  will  be  the  analysis, 
in  Part  IV,  of  a  simply  supported  spline  in  tension  from  a 
"Strength  of  Materials"  viewpoint. 
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II.  SIMPLY  SUPPORTED  SPLINE 
2.1  Basic  Equations 

Let  the  ordinates  be  given  at 

respectively.  For  convenience  consider  the  unit  spacing 
x^...j^  -  Xj^  =  1  (i  =  0,  ...,n).  Let  y(x)  be  an  Interpolation  curve 
through  these  points  and  define  y^  and  y^'  as  the  first  and 
second  derivatives  of  y(x)  at  Xj^  for  1  =  0,l,..„,n+l.  Bie  Inter¬ 
polation  curve  associated  with  the  deflection  curve  of  a  simply 
supported  elastic  beam  (spline)  Is  made  up  of  a  sequence  of 
cubic  polynomials  fj^(x)  such  that 


y(x)  =  f^(x)  for  <  X  <  x^^^^  (i  =  0,  ...,n) 


yK+i>  =  ,  ^n^*^ 

*  "^^+1 


(1) 


with  the  following  conditions  on  value,  slope,  and  second 
derivative ; 


fl(Xi)  = 

=  y^^  1  =  0,  ...,n. 

(2) 

11m 

X  ->Xj^ 

fi.l(x) 

=  y^  1  =  1, .  ..,n-!-l. 

(3) 

lim 

X  ->x^ 

f{.l(x) 

=  ^^(xj^)  1  =  1,  ...,n. 

(4) 

11m 

X  ->Xj^ 

=  fi'(x.)  1  =  l,...,n. 

(5) 

The  cubic,  f^(x),  can  be  expressed  In  terms  of  the  ordinates 
and  slopes  at  each  of  Its  end  points,  x^^  and  : 
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fj^(x)  ■=  +  (x-x^)y^^.^  + 

{xi+i-x)(x-Xi)  {[y^  -  (yi+i’J'i)  J(xj^^.^-x)  (6) 

-  [y^+i  -  (yi+i-yi)](x-x^)} 
for  Xj^  <  X  <  Xj^^^  (i  =  0, . .  .,n) . 


The  form  of  fj^(x)  shovm  in  (6)  identically  satisfies  the  con¬ 
ditions  on  value  and  slope  expressed  in  (2),  (3)  and  (4);  however, 
it  does  not,  in  general,  satisfy  the  continuity  condition  on 
the  second  derivative  expressed  in  (5) .  In  order  to  expand  (5) 
v/e  first  differentiate  (6)  twice: 

fi'(x)  »  y^+j^  -  y^  +  3[2(yj^^j^-yj^)  -  ][xj^+;,^-2x+Xj^] 

for  Xj^  <  X  <  (1  ==  0,  ...,n). 

Prom  (7), 

f^(x^)  =  y|+3^-y^+3[2(yj^^.^-yj^)  -  (yi+i+yj))  (l  =  0,  ...,n)  (8) 


lim 

X  ->x^  ^  ^ 


yi-yj[_i-3[2(yi-yi_i)-(3’^^+yi_i)  J  (1=1, . .  .#n+i) . 

(9) 


Using  (8)  and  (9),  (5)  becomes 

y^.l  +  4y^  +  y^^j^  =  3(yj^+]^  -  y^.j)  (i  =  l,...,n).  (lO) 

Equation  (lO)  gives  n  equations  for  the  n+2  unknowns  y^  (l  =  0, 

...,nTl)  necessary  to  determine  y(x)  throughout  x^  <  x  _<  x^^^^ 

(cf.  (1)  and  (6)).  The  two  additional  required  equations  are 

usually  obtained  by  specifying  end  conditions  on  slope  or 

second  derivative  at  x^  and  x  . , . 

o  n+1 
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If  the  end  slopes  are  specified,  the  nxanber  of  unknowns 
is  immediately  reduced  to  n  and  these  will  be  determined  by  (lO) 
alone . 

If  the  second  derivative  y^  is  specified  instead  of  y^ 
then  the  equation 

or,  using  (7)  to  expand  (ll)  into  a  more  convenient  form, 

2y<5  +  y;{  =  -iyo  +  3(yi-yQ)  (12) 

is  added  to  eqs.  (lO)  to  provide  n+1  equations  for  the  n+1 
unknowns,  y^  (l  *=  0, ...,n).  Likewise,  the  equation 

y"+l  =  llm  f"(x),  (13) 

or,  using  (7)  to  expand  (13)* 

K  2y^+i  =  K+i  ^  3(yn+i  -  y^)*  (i4) 

is  considered  v^lth  eqs.  (lO)  when  y^^j^  is  specified  instead  of 
If  the  physical  spline  is  simply  supported  at  the  end¬ 
points,  the  internal  moment  at  those  points  is  zero.  In  this 
idealized  representation  of  the  physical  spline,  the  specifica¬ 
tion  of  y^  *  y^^j^  =  0  is  equivalent  to  having  zero  end  moments 
and  is  the  most  frequently  used  end  condition. 

2 .2  Conditions  for  an  Inflection  Point 

Assume  that  am  interpolation  curve,  y(x),  based  on  the 
above  cubic  polynomial  representation  of  a  spline,  has  been 
determined.  We  wish  to  determine  the  conditions  which  artt 
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assoclated  with  the  appearance  of  inflection  points.  Since  fj^(x) 
is  specified  by  its  endpoint  values  and  slopes  y^, 

the  einalysis  of  one  interval  will  yield  results  typical  of  all 
Intervals.  Consider  the  Interval  <  x  <  Since  f^'(x) 

(cf.  (7))  is  linear  there  can  be  an  inflection  point  if  and  only 
if  y^'  and  are  of  opposite  sign.  In  view  of  the  discussion 

in  Part  1,  we  must  study  the  second  differences  associated  with 
Xj^  and  in  order  to  determine,  assuming  y^  and  are  of 

opposite  sign,  if  such  an  inflection  point  is  extraneous .  Before 
doing  this  it  will  be  helpful  to  discuss  the  inflection  point 
conditions  in  terms  of  slopes  rather  than  second  derivatives. 

Let  s(x)  be  the  difference  between  the  slope  of  y(x)  and 
a  polygonal  function  connecting  the  given  points.  Then,  in  the 
internal  being  considered, 

8(x)  =  y«(x)  -  -  y^^)  <  x  <  (15) 

with  its  values  at  x^  and  being 

^  =  yi  -  ^^1+1  -  yi^ 

and 

°i+i  ~  yi+1  ■  ^yi+1  "  yi^' 

respectively.  Substituting  (16),  (17)  in  (8)  and  (9),  the  second 
derivatives  y^'  and  y^^^  can  be  expressed  in  terms  of  s^  and  : 

K  ^  “  ®1+1  "  ®i  ■  3(3i+i  Si)  (18) 

yi+1  ”  fj^(x)  =  Sj^^j  “  ®i  3^®i+i  ®i^ 

*“^*1+1 


(19) 
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Deflnlng 


(20) 


the  ratio  of  second  derlvates  has  the  form 


l+2Ci 


(21) 


this  ratio  is  positive  if  and  only  if  -2  <  Cj^  <  -  i*  If  and 
have  the  same  sign  (Cj^  >  O)  we  would  expect  an  inflection 
point.  For  s^  and  of  opposite  sign  an  inflection  point 

occurs  if  and  s^^^^  differ  in  magnitude  by  more  than  a  factor 
of  two;  an  inflection  point  in  this  case  is  not  usually  desirable. 
For  example,  consider  th'i  single  Interval  problem  (n  =  O)  shown 
in  Figure  2.  With  y^  =  0,  y^^  =  and  specified  and  slopes 
y<J  =  1,  y{  0,  we  have,  from  (l6),  (17)  and  (20),  s^  = 

5^  =  and  Cq  =  “3.  Tlie  cubic  polynomial  interpolation  curve, 
represented  by  the  solid  line,  shows  an  extraneous  inflection 
point  and  subsequent  reverse  in  curvature  near  the  right  end. 
Clearly,  a  curve  such  as  the  broken  curve  in  Figure  2  has  the 
desired  smoothness  which  the  solid  curve  lacks,  even  though  its 
second  derivative  must  have  a  larger  raaxlmiun  and  could  be  con¬ 
sidered  less  "smooth"  in  a  mathematical  sense.  Ihls  situation 
suggests  that  some  additional  condition  on  curvature  vfould  be 
useful  to  further  define,  mathematically,  the  desired  type  of 


smoothness . 
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III.  CONDITIONS  CN  CURVATURE. 


In  order  to  avoid  extraneous  inflection  points,  the 
discussion  concerning  Figure  1  leads  to  the  requirement  that  y^' 
have  the  same  sign  as  the  second  difference  at  which  will 
be  denoted  as  dj^.  Although  this  brief  statement  is  intuitively 
plausible,  it  leaves  two  definitions  and  some  Intervening  argu¬ 
ments  in  need  of  elaboration. 

Wie  first  problem  is  that  the  usual  definition  of  the 
second  difference  at  x^, 

^i  =  (y^^+i-y^)  -  (yi-yi_i)  =  yi+i  -  2yj^  +  yi_i  (i  =  i,...,n 

(22a) 


leaves  us  without  an  indicator  of  the  desired  sign  of  y|'  at  the 
endpoints  x^  and  If  y^'  is  specified  then  the  sign  is  ob¬ 
viously  the  desired  one,  however,  for  notational  convenience  in 
stating  the  curvature  conditions,  we  formally  define: 

do  =  y^,  if  y^  is  specified. 


and  similarly 


(22b) 


'^n+1  =  yn+1-  ^n+l  specified. 

For  the  case  of  a  specified  end  slope  vie  look  again  vo  Figure  1. 
If  the  specified  slope  at  1  is  that  of  the  dashed  line  a,  we 
would  expect  a  positive  second  derivative  at  Ij  if  the  slope  is 
that  of  b  then  a  negative  second  derivative  is  desirable.  A 
;onvenient  Indicator  which  duplicates  these  results  is  a  second 
difference  which  uses  the  specified  slope  as  the  leading  first 
difference;  l.e. 
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d©  =  “  ^0'  specified, 

and  similarly  (22c) 

Vl  =  ^n+l  -  ^yn+1  ■  K+1  specified. 

Definition.  Assume  that  y(x)  is  an  InteiTpolation  curve 
of  class  ^  i5CQ»*n+l^  through  the  given  points  at  (l  »  0, 
...,n+l),  satisfying  conditions  on  slope  or  second  derivative 
at  and  Assume  further  that  y(x)  has  at  most  one  in¬ 
flection  point  in  the  interval  Xj^  <  x  <  (i  ■=  0,  ...,n). 

Then,  such  an  inflection  point,  if  it  exists,  will  be  extraneous 
if  dj^  and  have  the  same  sign. 

Theorem  1.  Assume  y(x)  as  in  the  above  definition. 

Then  the  condition  that  yj|  and  dj^  (l  =  0,  ..,,n+l)  have  the  same 
sign  is  necessary  and  sufficient  for  y(x)  to  have  no  extraneous 
Inflection  points . 

Proof .  The  condition  is  clearly  sufficient  since  for 
every  interval  such  that  d^^  and  have  the  same  sign,  y^'  and 

^i+1  have  the  same  sign  and  there  vjould  be  no  inflection 

point.  Proof  of  necessity  is  more  involved  and  will  be  demon¬ 
strated  with  the  aid  of  Figure  3  and  the  notation  that,  say, 
(+,-)j^  will  mean  is  positive  and  is  negative.  Assume  that 
there  is  some  point  i  such  that  y^'  and  d^^  have  the  same  sign, 
say  (+,  +  )j^,  and  and  dj^^^  have  opposite  signs.  We  have  two 
choices,  ^the  former  implies  an  extraneous 

Inflection  point  so  we  choose  the  latter,  iecause  of  the  signs 
of  dj^  and  it  is  clear  from  Pigxire  3  that  y^'+2 
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negatlve  since  only  one  inflection  point  in  an  interval  is  per¬ 
mitted;  then  only  comLination  that  doesn't  make 

the  inflection  point  extraneous .  Repeating  the  argvunent  for 
i+3>  1+^,...  we  find  that  the  alternation  of  signs  must  continue 
if  extraneous  inflection  points  are  to  be  avoided.  In  this  case, 
however,  y(x)  cannot  satisfy  the  final  boundary  condition.  For 
instance,  if  1+4  is  the  last  point,  the  dashed  line  would  be  a 
typical  specified  slope  associated  with  the  positive  a 

value  impossible  for  y^+jj  to  achieve.  If  hsd  been  specified 

the  contradiction  is  more  direct  since  d^^^;  the  same  sign  by 
definition.  The  initial  assumption  that  there  is  at  least  one 
point  1  such  that  dj^  and  y^  have  the  same  sign,  if  there  are  no 
extraneous  inflection  points,  can  be  shown  by  an  argument  similar 
to  the  one  Just  completed.  If  one  assumes  d^  and  y^  have  oppo¬ 
site  signs  at  all  points  then  the  boundary  condition  at  either 
the  first  or  last  point  cannot  be  satisfied. 

The  cubic  interpolation  curve  may  satisfy  the  above 
conditions  on  the  second  derivative  (i.e.,  y^'  and  d^^  have  the 
same  sign),  but  if  it  does  not,  nothing  can  be  changed  to  correct 
it.  In  order  to  find  a  spline  curve  which  has  a  parameter  that 
can  be  varlea  to  satisfy  the  curvature  conditions,  we  are  led, 
with  some  physical  justification  as  vrell,  to  consider  a  simply 
supported  spline  in  tension. 
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IV.  SIMPLY  SUPPORTED  SPLINE  IN  TEi'ISION 
4.1  Basic  Equations 

Por  the  spline  in  tension  we  have,  corresponding  to  (6), 
the  form  [8]; 


f^(x)  =  +  (x-x^)yj^^^  + 

+ -^[y^  ■^’1^1+1  ■  (n+i)(yi+i-yi)  H 

n  -1 


sinh  p(x-x^)-(x-Xj^)sinh  p 
sinh  p  -  p 


] 


sinh  p(xj^^j^-x)-(x^^, -x)sinh  p 


1+1 


], 


for  Xj^  <  X  <  Xj^+2^ 


where 


and 


sinh  p  -  p 

(l  =  0, ...,n), 

^  p  cosh  p  -  sinh  p 
sinh  p  -  p 

P  oc  ^/tension 


(23) 


(24) 


Although  not  immediately  obvious,  (23)  reduces  to  (6)  for  p  =  0. 
Figure  4  shows  a  graph  of  rj  as  a  function  of  p.  Por  p  large,  t) 
approaches  p-1  asymptotically. 

Co/cresponding  to  (?),  the  second  derivative  is 

,  f  sinh  o(x-x. ) 

fi'(x)  =  ^nyl+i  ■  (n+i)(yi+i-yi))-  ^ 


-  [nyi  +  y;[+i  -  (n+i)(yi+i-yi) ) 


sinh  p  -  p 

p‘^  sinh  p(xj^^.^-x)  I  ^ 


for  Xj^  <  X  < 


sinh  p  -  p 
(1  =  0, ...,n). 


J 


(25) 
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Prom  (25), 


^[fiyi+yi+i-(n+i)(yi+i-yi)  ] 

T)  -1 


2 

p  slnh  p 
sinh  p  -  p 


(26) 


and 


lira  f!'  ,(x) 


4-^fyi_i''’TyjL~(n+i)(yi-yi.i)  I 


2 

p  sinh  p 
sinh  p  -  p 

(27) 


As  In  the  cubic  polynomial  representation  of  the  spline,  the 
form  of  fj^(x),  (23) i  assures  that  the  conditions  on  value  and 
slope,  expressed  In  (2),  (3)  and  (4),  are  Identically  satisfied. 
Using  (26)  and  (2?)/  the  continuity  condition  on  the  second 
derivative,  (5)^  becomes 


y^.l  +  2r)y^  -i-  y^^j^  =  (n+i)  (yi+i-y^.i)  (1  =  1,-.  .^n)  (28) 


Note  that  since  r|  equals  2  In  the  limit  as  p  approaches  zero, 

(28)  reduces  to  (lO)  as  p  0. 

Equation  (28)  gives  n  equations  for  the  n+2  unknowns 

y^,  (1  =  0,1, . .  .,n+l) ,  If  the  end  slopes  y^  and  are 

Ep3. If Lcd,  there  are  only  n  unloiowns  and  equations  (28)  are 

sufficient.  If  the  second  derivatives  at  the  ends,  y"  and 

o 

y^^^  are  specified,  then  the  equations  (cf.  (12),  (l4)) 


qy'  +  y-l  =  -y 


sinh  P  -  p  /  2  , \  /  •  /  \ 

0  — 2— (n  -1)  +  (n+x^'y^-:.' ) 
°  p  sinh  p  ° 


y'+qy/ni  =  4m  (n+i)(yn+T.-yn), 


(29) 


P  sinh  p 

are  combine(i  v;ith  (28)  to  provide  n+2  equations  for  the  n+2 


unknowns . 
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4.2.  Conditions  for  an  Inflection  Point. 

Combining  (25)  with  (26)  and  (27),  the  second  derivative 
can  be  expressed  In  terms  of  its  values  at  the  Interval  end¬ 
points  : 


.  sinh  p(x. ,,-x)  sinh  p(x-x. ) 

f "(x)  =  v"  ^  - _ -  +  v"  _ — _ _ 


sinh  p 


<  X  < 


sinh  p 

(l  =  0^...jn) 


(30) 


Inspection  of  this  form  shows  that  although  f^(x)  is  not  a 
linear  function  as  in  the  cubic  polynomial  case.  It  still  has 
the  property  that  it  will  not  cross  zero  unless  the  end  values, 
y^'  and  have  opposite  signs.  For  instance,  if  y^'  and 

^1+1  positive  then,  going  from  x^^  to  the  first 

term  in  (30)  decreases  monatonically  frcm  yj|  to  zero  while  the 
second  tern  increases  monatonically  from  zero  to 
both  terms  are  positive  In  the  Interval  and  there  can  be  no 
inflection  point.  A  similar  argument  shovfs  that  if  yj|  and  y^'+j^ 
are  of  opposite  sign  then  there  is  one  and  only  one  Inflection 
point.  Thus,  we  have  shown  that  the  interpolation  curve,  y(x) 
(cf.  (l)  -  (5)),  using  the  hyperbolic  spline  representation 
f^(x)  shown  in  (23) j  satisfies  the  assumptions  on  y(x)  specified 
in  the  Definition  and  Theorem  1  discussed  in  Part  III. 

The  inflection  point  conditions  in  terms  of  s^^  are 


developed  in  same  manner  as  that  discussed  In  Part  2.2  for 
the  cubic  polynomial  representation  (see  (15)  -  (2l)).  Bie 
final  results  are  that  there  will  not  be  an  inflection  point 
in  the  Interval  Xj^  <  x  <  if  and  only  if  -f]  <Cj^  <  -  I/4 
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and  that  for  and  of  opposite  sign,  an  extraneous  in¬ 
flection  point  v.'ill  not  occur  if  s^  and  differ  in  magni¬ 

tude  by  less  than  a  factor  of  t)  .  At  this  point  one  may  be 
led  to  the  conclusion  that  by  increasing  tension  sufficiently, 
will  accommodate  any  specified  ratio  C  and  therefore  all 
extraneous  inflection  points  can  be  removed.  Ihe  problem  is 
that  3^  is  a  function  of  »i  and  thus  as  changes,  so  Goes 
Since  the  hyperoollc  spline  representation  approaches  a  poly¬ 
gonal  function  as  fj  becomes  large,  the  values  of  K,  ^  approach 
finite  limits,  which  places  the  conclusion  on  firmer  ground. 

It  is  more  useful,  however,  to  establish  this  conclusion  by 
completely  determining  the  effect  of  tension  on  slope  and  curva¬ 
ture  throughout  all  the  intervals.  Clearly,  an  explicit  solu¬ 
tion  for  the  unlcnowns  y^  (l  =  0,  ,.,n+l)  is  desirable. 


4.3  An  Explicit  Solution. 

Consider  first  the  case  of  given  y^  and  ^n+l*  Equa¬ 
tions  (28)  can  be  written  in  matrix  form, 

N  Y'  =  C  (31) 


with 


2fl  1  0  . 

1  2n  1. 

0  1  2t]' 

«  « 

•  < 

•  % 

0 

0 

0  0  .  .  . 


,.000 

0 

«  • 

•  ^  • 

/2ti  *1  0 

1  2ri  1 

.  0  1  2ti 


(32) 
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I  V» 


yi  = 


y; 


n-] 


y  * 


(33)  i(T'  l)(y2-yo^  -  |(34) 

(n+i)(y3-yi) 

(n+i)(y4“y2) 


c  = 


(ri+i)(yn.i-yn_3) 

(n+i)(yn-yn_2) 

(l+l)(yr,+i-yn-i)  " 

since  n  ^  2  for  all  p  >  0,  the  diagonal  eler.icnt  in  each  row 
of  N  is  larger  than  the  sum  of  the  absolute  values  of  the  re¬ 
maining  elements  in  that  row;  thus,  by  a  well  known  theorem 
on  determinants  [9],  N  has  an  inverse.  In  fact,  the  Inverse 
of  N  (say,  order  ir.)  can  be  written  explicitly  [10,11]: 


r 


N'^=  ^ 
m 


p 

P  T 

-p  p  0 

+P  P  0  . 

0 

m-1 

0  m-2 

0  m-3 

p 

P  ^ 

+P,P_  « 

-P,P  ,  . 

0 

m-2 

±  m-t: 

j.  m-j 

p 

P  0 

-P,P  , 

+P^P  0  . 

0 

m-3 

1  ra-3 

2  m-3 

kin+l, 


o  o 


im 


>m+l 


1*^0 

P  ? 
^2^0 


nr^l 


where 


^o  =  1 

=  2  n 

Pg  =  2r)P,  - 

P3  ~  2  t]  Pg  -  Pj^ 

Pk  =  2  n  P,  -  Po 


(35) 


^m  =  2ti  P„_i-  P„, 


m-1  *ra-2 


(36) 
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or,  since  is  the  Tschebyscheff  polynomial  of  the  second  kind^ 
U^(r;),  the  direct  representation  is 


F.  =  2"^  (  T]  -  cos 

1:-1 


•j:^)  (i  - 


"o  = 


Defining  a  matrix  R  as  the  adjoint  of  N, 


(37) 


R  =  P  N' 


(38) 


equation  (35)  can  be  written  In  Indlclal  form: 

\j  -  fi-iVj  1  i  J  ‘-J  - 1 . ” 


(38') 


.11  iJ 


Ijj  l^»..^m 


Also  we  note  symmetry  about  the  non-prlnclpal  diagonal: 

^‘m+l-j  m+l-l  "  hj  ^  l,...,m 

In  the  above  case  of  given  end  slopes,  m  =  n.  In  the 
case  of  given  y^  and  y^+j  given,  N  takes  on  a  slightly  changed 
form  since  equations  (29)  must  be  Included: 


0  .  . 


1  2n  1 
0  1.  2n  '  . 


0  0  0 


':39) 


0  0  . 


2t)  '1  0 

'1  2f)  1 

0  1  n 
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with  order  m  =  n+2.  The  inverse  of  this  second  form  of  N,  is 
identical  to  the  form  shown  in  (35)  if  is  redefined  as  [12,11] 

^o  =  1 
Pi  =  n 


P2  =  2nPi  -  Po 

Pj  =  2T1P2  -  p^ 


(40) 


Pm  1=  2flP_  o  -  P„  o 

in-2  in-»j 

P_  =  T  Pm  1  ~  Pm  O 

m  '  m-i  ni-2 

Except  for  P^j^,  the  P^  (i  =  0,  ...,m-l)  in  (40)  are  Tschebyscheff 
polynomials  of  the  first  kind,  Tj^(r)),  and  have  adirect  represent¬ 


ation 


"0  =  ' 

P.  =  2^'^  n  (fi-cos%ix)  (i  =  1,  ...,m-l) 


With  y^  and  y^+i  specified,  the  column  vectors  Y’  and  C 
(cf.  (31))  that  correspond  to  N  (cf.  (39))  are; 


r^o  ] 


(42) 


y"  = 


^n+l 


(ti+i)(yi-yo)  -  (n^-i)  ^ 

S  j 

(n+i)(y2-yo) 


p  sinh  p 


(n+i)(y3-yi) 

(n+l)(y^l-yn.i) 

(n+i)(yn+i-yn)  +  (r-i)"z' ;  .  y" 

n-rj.  n  p^’sinh  n 


(43) 
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Because  of  the  similarity  In  the  forms  of  the  solution,  the  case 
of  given  and  y^^j^  will  not  be  discussed  further. 

Theorem  2.  Let  an  Interpolation  curve  y(x)  be  of  the 
r--rm  specified  In  eqs.  (l)  -  (5)  with  fj^(x)  (l  =  0,  ...,n)  being 
the  hyperbolic  spline  representation  given  in  (23) .  Assume  that 
the  end  slopes  and  are  specified  and  that  dj^  /  0  (l  =  1, 

...,n),  (For  the  case  of  dj^  =  0,  l.e.,  three  given  points  In  a 
l:',ne,  s<=e  Part  4.4,)  Then,  the  condition  expressed  in  (5)  and 
equivalently  In  (28)  and  ln(3l)has  the  solution 

Y'  =  N“^C  (44) 

where  N  ^is  given  by  (35).  Thus,  y(x)  Is  completely  determined. 

Further,  there  exists  a  value  of  t)  ^  2  such  that  the 
curvature  conditions  expressed  In  Itieorem  1  are  satisfied,  l.e. 

y;;  d^  >  0  (1  =  0,...,n+l).  (45) 

Proof.  First  we  will  show  that  (35)  is  the  ccrrect  inverse. 
Using  (38),  this  Is  equivalent  to  proving. 

NR  =  PI  (46) 

B  =  NR  (47) 

then  using  (32)  and  (38')  the  diagonal  terms  can  be  expressed: 

(1  7  l.m) 

®11  ^■^1-2'^m-l  ^^^l-l^m-l"^l-l^m-l-l^ 

(48a) 

=  p,  ,(2np„  ,  -  p„  ,  ,)  -  p,  iP„  .. 

1-1'  m-i-1  1-1  m-1 
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Using  (36)  repeatedly: 


®11  ~  ’^1-1  ^m-l+l  "  ^1-2  ^m-1 


(2t)P^_2  -  ^i-3^^ra-i+l  ■  ^l-2^ra-i 
=  Pl_2(2nPn,_i_+i  - 


=  P  P  -  P  P 

^l-2^m-l+2  ^1-3  m-1+1 


=  P  P  -  P  P 

m-1  o  m-2 


Also  B. 


"  ^m-2  “ 

"  "  ^o^m-2  * 

“  ^m-1  ■  ’^m-2 


(48b) 

(48c) 


Since  N  is  symmetric  we  need  show  only  the  upper  matrix  terms 
( J  >  1 )  as  zero : 

(j  >  1,  1  7  1) 


(J  >  1,  1  =  1) 


=  (-1)^  ^2n)]  =  0 


(48d) 


(48e) 


Thus,  eqs.  (48)  verify  (46). 
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Nexti  consider  the  existence  of  a  value  of  r)  >  2  such 
that  there  are  no  extraneous  Inflection  points .  Ihe  second 
derivative  at  1,  can  be  expressed  in  two  equivalent  (cf.  (5)) 
forms,  (26)  and  (27).  By  adding  these  two  forms  we  obtain  the 
more  symmetric  expression: 


2yi  =  7-I— -  -  y^+i  +  (n+i)(yi+i  -  sy^  + 

(q  -1)  s.'.nh  p-p 

(l=l,...,n)  (49) 

or,  using  (22) 

y;'  - 5— -  .  -  y.  +  (r)+l)d  ]  (1  =  l,...,n)  (49') 

^  2{rr-l)  slnh  p-p  ^  ^  ^ 

Substituting  (49')  into  (45)  and  omitting  the  positive  constants, 

v/e  must  show  that 


[y^.l  -  y^+i  +  (ri^l)djdj  >0  (i  =  0,...,n+l)  (50) 

for  some  tj  >  2 . 

Assume  >  0  and  consider  the  1^^  inequality  of  (50) 

(i  /  0,n+l): 

(n+i)d^  >  y^+i  -  y^.i-  (5i) 

Using  the  solution  for  y^  (cf .  (44))  where  N“^,  (35)>  is  expressed 
in  terms  of  the  elements  cf  its  adjoint  (38),  the  inequality  (51) 
becomes : 


(n+i)d,  Pj,  >  2 

*  "  j=i 


(p 


1+1  J 


From  (36)  we  see  that  is  an  i^*^ 
(36)  and  (38'),  (1,J  =  1, ...,n) 


(52) 


order  polynomial  in  t|  ;  from 
is  a  polynomial  of  order 
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n-1  or  less;  frcn:  (34),  is  linear  in  t]  .  Further,,  it  :;an 
be  shown  from  (36)  that 

0<  Po<  Pj  <P2  <  ...  <p„,  ni2.  (53) 

Thus,  we  find  that  the  R....S.  of  (52)  is  a  polynomial  in  t]  of 
order  n  while  the  L.H.S.  of  (52)  is  a  polynomial  in  t)  of  order 
n+1  and  is  positive  for  tj  >  2.  Thus,  there  will  always  be  a  value 
of  n  >  2  which  satisfies  the  Inequalities  (52) .  Such  a  value 
v.'ill  be  the  largest  root  of  a  polynomial  of  order  n+1  and  can  be 
easily  found  ntunerlcally  if  not  analytically.  Other  cases  in 
the  proof  follov;  essentially  this  Scime  argument. 

Having  determined  the  minimum  value  of  which  satisfies 
the  curvature  conditions,  the  unknown  slopes  can  be  determined 
from  (44);  (23)  is  then  used  to  perform  the  Interpolation  in  the 
desired  regions . 

An  example  of  a  "tensioned"  spline  curve  is  shown  as  a 
broken  line  in  Figure  5.  'Rie  solid  line  shows  a  normal  spline 
curve  through  the  3  given  points  0;  1.0,  2.0.  Ihe  two  intervals 
shown  are  actually  part  of  a  7  interval  problem;  the  vertical 
scale  is  exaggerated  5  times  for  clarity.  No^e  that  two  extrane¬ 
ous  inflection  points  are  removed  by  the  use  of  the  spline  in 
tension . 

4.4  Some  Practical  Considerations 

In  the  case  where  three  consecutive  points  are  in  a 
straight  line,  i.e.  dj^  =  0,  the  conditions  (45)  cannot  be  satis¬ 
fied  for  finite  q  .  An  infinite  value  of  tj  will  Tuake  the 
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interpolation  curve  simply  a  polygonal  function  -  which  is  harcly 
satisfactory.  Practical  considerations  suggest  prescribing  a 
straight  line  from  to  ^nd  considering  two  problems: 

a  curve  from  to  ^  with  specified  and  a  curve  from 

^i+1  ^n+1  ^1+1  specified.  Where  d^^  is  very  close  tc  zero, 

a  circular  segment  prescribed  through  the  given  points  at 
Yl,  and  the  use  of  the  above  scheme  vrould  probably  be  pref¬ 

erable  to  a  very  large  value  of  tj .  In  general,  specification 
of  a  discontinuity  in  curvature  or  slope  at  an  intermediate  given 
point,  or  an  interval  of  specified  curvature,  can  be  handled  by 
considering  two  spline  curves,  one  on  either  side  of  the  specified 
discontinuity  or  region. 

If  a  predominance  of  problems  deal  with  specified  end 
curvatures,  it  is  useful  to  derive  the  solution  in  terms  of  second 
derivatives  at  the  given  points  Instead  of  the  slopes  since  it 
avoids  the  explicit  use  of  the  two  additional  equations  (29) . 

In  such  a  change  of  variables,  (4)  rather  than  (5)  are  the  con¬ 
ditions  which  are  not  identically  satisfied,  but  the  resulting 
matrix,  N  (cf.  32)),  is  unchanged. 
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